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Abstract: The characterization of chaotic Hamiltonian systems in terms of the curvature 
associated with a Riemannian metric tensor in the structure of the Hamiltonian can be 
extended to a wide class of potential models of standard form through definition of a 
conformal metric. The geodesic equations reproduce the Hamilton equations of the orig- 
inal potential model when a transition is made to an associated manifold for which the 
geodesies coincide with the orbits of the Hamiltonian potential model. We therefore find a 
direct geometrical description of the time development of a Hamiltonian potential model. 
The second covariant derivative of the geodesic deviation in this associated manifold gener- 
ates a dynamical curvature, resulting in (energy dependent) criteria for unstable behavior 
different from the usual Lyapunov criteria. We discuss some examples of unstable Hamil- 
tonian systems in two dimensions giving, as a particular illustation, detailed results for a 
potential obtained from a fifth order expansion of a Toda lattice Hamiltonian. 
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A Hamiltonian system of the form (we use the summation convention) 

B = j^WkV, (1) 

where gij is a function of the coordinates alone, has unstable orbits if the curvature asso- 
ciated with the metric is negative. One can easily see that the orbits described by the 
Hamilton equations for (1) coincide with the geodesies on a Riemannian space associated 
with the metric g^, > i.e., it follows directly from the Hamilton equations associated with 
(1) that (using (12) and the time derivative of (10)) 

X£ = x m x n , (2) 

where the connection form VJ 

■pmn 

1 e 

and g^ is the inverse of g^. 



is given by 

1 ,dg km dg kn dg nm , 
2 gek l g x q x q x 
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The second covariant derivative of the geodesic deviation depends on the curvature 2 ' 3 



1 ~ dxi dx k m 1 m 1 ' 1 } 

i.e., for £i = x\ — Xi on closely neighboring trajectories at t, 

J^T = R J / k XjX k &, (5) 

where D/Dt is the covariant derivative along the line Xj(t). The sign of the scalar con- 
traction of (4) then gives information on the stability of the orbits. 3 

In this Letter, we point out that this formulation of dynamic stability has application 
to a much wider range of Hamiltonian models; in fact, every potential model Hamiltonian 
of the form 

H=£ + V(x), (6) 

where V is a function of space variables alone, can be put into the form (1), where the 
metric tensor is of conformal form. 4 We obtain in this way a direct geometrical description 
of the time development for a Hamiltonian potential model. 

Casetti, Pettini and collaborators 5 , for example, have studied the application of both 
the Jacobi and Eisenhardt metrics in their analyses of the geometry of Hamiltonian chaos. 
The Jacobi metric 1 (of the form (E — V)5ij) leads to geodesic equations parametrized by 
the invariant distance associated with this metric on the manifold, in this case, the kinetic 
energy, thus corresponding to the Hamilton action. Transformation to parametrization by 
the time t leads to the second order Newton law 5 in the form (14) below, for which the 
geometrical structure is no longer evident. 

The Eisenhardt metric, leading to geodesic motion in t, involves the addition of an 
extra dimension. As noted by Caini et a/ 5 , this metric leads to the tangent dynamics 
commonly used to measure Lyapunov exponents in standard Hamiltonian systems. The 
method that we use, associated with a curvature that is explicitly energy dependent, 
appears to be a more sensitive diagnostic than the computation of exponents of a locally 
linearized system. 

The formulation of Hamiltonian dynamics of the type of Eq.(6) in the form (1) is 
carried out by requiring that (6) be equivalent to (1). For a metric of conformal form 

9ij = (7) 

on the hypersurface defined by H = E = constant, the requirement of equivalence implies 
that 

Substituting this result in the geodesic equations (2), one obtains an equation that 
does not coincide in form with the Hamilton equations obtained from (6). 
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To see that the Hamilton equations obtained from (1) can, however, be put into 
correpondence with those obtained from the Hamiltonian of the potential model (6), we 
first note, from the Hamilton equations for (1), that 

dH 1 

We then use the geometrical property that X{ is a first rank tensor (as is p l ), under local 
diffeomorphisms that preserve the constraint that H be constant, to define the velocity 
field 

^5^ = ^', (10) 

coinciding formally with one of the Hamilton equations implied by (6) . From this definition, 
we recognize that we are dealing with two manifolds, each characterized, as we shall see, 
by a different connection form, but related by 

dx j = g Jl dx z (11) 

on a common tangent space at each point (for which g^ is nonsingular). 

To complete our correspondence with the dynamics induced by (6), consider the 
Hamilton equation for p l , 

With the form (7) for g % \ we obtain in the particular coordinate system in which (6) is 
defined, 

p= -e=vW,- (13) 

Considering (11) as a change of variables, (13) becomes 

the second Hamilton equation in the usual form, where V is considered a function of the 
{x }, now considered as independent variables. 

As a coordinate space, the {x e }, which we shall call the Hamilton manifold, is not 
uniquely defined in terms of the original manifold {xi}, which we shall call the Gutzwiller 
manifold, since (11) is not an exact differential. As we have remarked, we shall be working 
with two manifolds (characterized by the connection forms (3) and (21)). It is the local 
relation (11) which induces, from the geometry of the Gutzwiller manifold, a corresponding 
geometry on the Hamilton manifold. We shall discuss applications and interpretation of 
the physics of the Gutzwiller manifold elsewhere, but turn now to a further examination 
of the consequences of the relations (10) and (111). 
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The geodesic equation (2) can be transformed directly from an equation for Xj to an 
equation for x\ the motion defined in the Hamilton manifold. From (10) it follows that 



xi = gijX J + -^-x n x J 

1 ( dg km dg kn dg nm , 
= ~2 9ek ^~dF~ + dx~ " ~d^! XmXr 

& vjjn v<L m vjyk 



(15) 



Now, using the identity 

dgij dg km 

" -9£k^—9mj, (16) 



re 



it follows that, with the symmetry of x n x m , 

— r x n x> = -^9ik{-^— + -w—)x n x m . (17) 

Thus, the term on the left side of (15) containing the derivative of gej cancels the first two 
terms of the connection form; multiplying the result by the inverse of gij, and applying 
the identity (16) to lower the indices of g nm in the remaining term on the right side of 
(15), one obtains 

x l = -M l mn x m x n , (18) 

where 

M^V'^T- (19) 
mn 2 dx k 

Eq. (18) has the form of a geodesic equation, with a truncated connection form. In fact, it 
can be shown (a full proof will be given elsewhere) that the form (19) is indeed a connection 
form, transforming as 

, rlP dx' e dx p dx q , dx' 1 d 2 x r 
M = M + 



mn 9x r dx' m dx' n pq dx r dx' m dx ,n ' 

consistent with the covariance of (18) under local diffeomorphisms of the Hamilton mani- 
fold. 

Substituting (7) and (8) into (18) and (19), the Kronecker deltas identify the indices 
of x m and x n ; the resulting square of the velocity cancels a factor of (E — V) -1 , leaving the 
Hamilton-Newton law (14). Eq. (18) is therefore a covariant form of the Hamilton- Newton 
law, exhibiting what can be considered an underlying geometry of standard Hamiltonian 
motion. 

The geometrical structure of the Hamilton manifold can be understood as follows. 
Let us write the covariant derivative for a (rank one) covariant tensor on the Gutzwiller 
manifold (defined as transforming in the same way as d/dx m ), using the full connection 
form (3), 

Q Am 

A m ' q = -T™ q A k . (20) 

8x q K 
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Lowering the index q with gi q , we obtain the covariant derivative in the Hamilton manifold, 
with connection form (with the help of (16)) 



F m _ „ r mq 1 n mq ( ^9iq dg k q &g M y 

r mk = 9eq r k =-g l^-^-^}- (2i) 

This induced connection form, in the formula for curvature, would give a curvature cor- 
responding the the Hamilton manifold. However, it is antisymmetric in its lower indices 
(£, k) (torsion). Taken along a line parametrized by t, corresponding to geodesic motion, 
the antisymmetric terms cancel, leaving precisely the symmetric connection form (19) * . 
A complete discussion of the tensors on the Gutzwiller and the Hamilton manifolds will be 
given elsewhere. We note here, however, that the curvature associated with the geodesic 
deviation in the Hamilton manifold, as we shall see below, is not the same as the intrinsic 
curvature of that manifold, determined by T^ fc , but rather, due to the presence of torsion, 
a special curvature form associated with the geodesies themselves. 

Since the coefficients M^ n constitute a connection form, they can be used to construct 
a covariant derivative. It is this covariant derivative which must be used to compute the 
rate of transport of the geodesic deviation ^ = x fi — x* along the (approximately common) 
motion of neighboring orbits in the Hamilton manifold, since it follows the geometrical 
structure of the geodesies. 

The second order geodesic deviation equations ** 

f) ft/ft- 

e = -2Mi n x m C - -^iVf, (22) 
obtained from (18), can be factorized in terms of this covariant derivative, 

efn = % + Mir. (23) 



One obtains 



= R M \mnX q X n r, (24) 



where the index M refers to the connection (19), and what we shall call the dynamical 
curvature is given by 



e dM* m dM* n 

?™ dx n 8x n (25) 



+ M k qm M" nk - M k qn Ml k . 



* Note that since (19) and (21) are not directly derived from g^-, they are not metric 
compatible connections. However, performing parallel transport on the local flat tangent 
space of the Gutzwiller manifold, the resulting connection, after raising the tensor index 
to reach the Hamilton manifold, results in exactly the "truncated" connection (19). 

** Substituting the conformal metric (7) into (22), and taking into account the constraint 
that both trajectories x' e and x e have the same energy E, one sees that (22) becomes the 
orbit deviation equation based on (14). 
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This expression, as remarked above, is not the curvature of the Hamilton manifold (given by 
this formula with in place of M^ m ), but a dynamical curvature which is appropriate 

for geodesic motion. 

We give in the following a general formula for the geodesic deviation in the Hamilton 
manifold in two dimensions, and then show results of computer simulation for Poincare 
plots showing a correspondence with the prediction of instability from the geodesic devia- 
tion. 

With the conformal metric in noncovariant form (7), (8), the dynamical curvature 
(25) can be written in terms of derivatives of the potential V, and the geodesic deviation 
equation (24) becomes 



- - (26) 



where the matrix V is given by 



, 3 dV dV ; 1 d 2 V 1 



M 2 v 2 dx t dx i Mdatdx* 

and 



v 1 v 3 

P x 3 = 8"-— ^, (28) 



with v l = x l , defining a projection into a direction orthogonal to v l . 
We then find for the component orthogonal to the motion 

Dll ^\ 2 = - [Ai cos 2 <P + A 2 sin 2 0] ( V± • (29) 

where Ai and A 2 are eigenvalues of the matrix V, and is the angle between and the 
eigenvector for Ai. 

Instability should occur if at least one of the eigenvalues of V is negative, in terms of 
the second covariant derivatives of the transverse component of the geodesic deviation. 

One may easily verify that the oscillator potential is predicted to be stable. Our 
criteria imply that the Duffing oscillator (without perturbation, not a chaotic system) 
clearly indicates instability in a neighborhood of the unstable fixed point. The potentials 
discussed by Oloumi and Teychenne 6 also demonstrate the effectiveness of our procedure; 
our results in these cases are in agreement with theirs. The relation (29) provides a clear 
indication of the local regions of instability giving rise to chaotic motion in the Henon- 
Heiles model (this result will be discussed in detail elsewhere). 

We take for a simple illustration here a slight modification of the fifth order expansion 
of a two body Toda lattice Hamiltonian (for which the fourth order expansion coincides 
with the Henon-Heiles model) 

V(x, y) = \(x 2 + y 2 ) + x 2 y - ^y 3 + + \y\ (30) 

This provides a new Hamilton chaotic system for which our criterion gives a clear local 
signal for the presence of instability. Fig.l shows that the region of negative eigenvalues 
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does not penetrate the physically accessible region for E = 1/6; fig. 2 shows a Poincare 
plot in the y,p y plane for this case, indicating completely regular orbits. In fig. 3, the 
distribution of negative eigenvalues for E = 3 is shown to penetrate deeply into the physical 
region, and fig. 4 shows the corresponding Poincare plot displaying a high degree of chaotic 
behavior. The criterion for instability we have given depends sensitively on the energy of 
the system. The critical energy for which the negative eigenvalues begin to penetrate the 
physically accessible region, in this example, is E = 1/5. 

The condition implied by the geodesic deviation equation (26), in terms of covariant 
derivatives, in which the orbits are viewed geometrically as geodesic motion, is a new 
condition for instability, based on the underlying geometry, for a Hamiltonian system of 
the form (6). This geometrical picture of Hamiltonian dynamics provides, moreover, new 
insight into the structure of the unstable and chaotic behavior of Hamiltonian dynamical 
systems. 

We wish to thank S. Shnider, A. Belenkiy, P. Leifer, I. Aharonovitch and Avi Gershon 
for helpful discussions. 
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FIGURE CAPTIONS 



Figure 1. 

The dark area shows the region of negative eigenvalues for the matrix V. The light area 
corresponds to physically allowable motion for E = 1/6. The region of negative eigenvalues 
does not penetrate the physically accessible region in this case. 

Figure 2. 

A Poincare plot in the (y,p y ) plane for E = 1/6, indicating regular motion. 
Figure 3. 

The dark area of negative eigenvalues for the matrix V is seen to penetrate deeply 
into the light region of physically allowable motion for E = 3. 

Figure 4. 

A Poincare plot in the (y,p y ) plane for E = 3, indicating strongly chaotic behavior. 
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